Information rate of a waveguide 
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We calculate the communication capacity of a broadband electromagnetic waveguide as a function 
of its spatial dimensions and input power. We analyze the two cases in which either all the available 
modes or only a single directional mode are employed. The results are compared with those for the 
free space bosonic channel. 
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In the analysis of electromagnetic communication 
channels using quantum information, emphasis has been 
placed on free space communication protocols, in which 
waves propagate unconstrained: a good summa of all the 
obtained results can be found in Here we will fo- 

cus on constrained communication lines, such as optical 
fibers or radio waveguides, in the lossless limit. Although 
the spatial mode structure for free space propagation be- 
tween a pair of apertures has long been understood Q , its 
near-field modes only approximate the lossless behavior 
of ideal waveguides. In contrast, since the spatial proper- 
ties of the waveguide modes are always well-defined, we 
will be able to derive the exact dependence of the infor- 
mation rate on the system parameters, e.g., the power 
P and the waveguide cross-sectional area A, obtainin 
results that closely resemble the ones described in 0, 
for the free space channel. Moreover, as will be discussed 
in detail, our derivation resolves some of the open issues 
connected with the optimization of the multimode com- 
munication protocols. 

We start by describing the waveguide communication 
channel in Secnjand calculate the rate in Sec.^ In par- 
ticular. Sees. Ill K\ and HlBl are devoted to the regimes of 
multiple modes and single directional mode respectively. 
The discussion and the comparison with prior results is 
given in Sect. IIIII 




FIG. 1: Description of the ideal metallic waveguide. The 
modes (TE or TM) with wave vector k propagate in the posi- 
tive longitudinal direction bouncing off the perfectly reflecting 
walls of the waveguide. 



I. THE CHANNEL 



Although guided-wave optical communications are 
normally carried out using dielectric waveguides, metal- 
lic waveguides provide a simpler mode structure for de- 
riving the broadband information rate. In the ideal, 
lossless case that we consider, such waveguides con- 
fine the electromagnetic field into a finite region of 
space by means of perfectly-reflecting boundaries. In 
this paper we will analyze in detail the rectangular 
cross-section case with transverse spatial dimensions Li 
and L2 described in Fig. ^ even though the proce- 
dure can be readily extended to other configurations. 
In a hollow waveguide the transverse electromagnetic 
modes (TEM) customarily used in free space commu- 
nications do not propagate. They are replaced by the 
transverse electric (TE) and the transverse magnetic 
(TM) modes in which only the electric or magnetic 
terms respectively possess null longitudinal components. 
These modes are characterized by the wave vectors 



ni, 712 = 1, 2, • 



TE 



Trni 7rn2 , 



ni,n2 = 0,1,2, ••• TM 
(ni — n2 — Q excluded) 

(1) 



where the discretization of the transverse components 
derives from the boundary conditions at the waveguide 
walls and where the longitudinal component is a pos- 
itive quantity because we are considering only modes 
propagating from the sender to the receiver. By intro- 
ducing the transmission time T (i.e. the time interval in 
which the sender operates the channel), the longitudinal 
parameter can be discretized using periodic bound- 
ary conditions. In particular, a mode with wave vector 
k bounces off the waveguide walls (see Fig. ^ so that it 
propagates across the transmission line with a longitudi- 
nal speed (group velocity) ccosO = ck^/\k\. This means 
that all the photons of that mode used in the transmission 
can be ideally enclosed in a box of longitudinal length 
cT cos 9: assuming periodic boundary conditions k^ can 
be discretized as 2imz/{cT cos 9) with 713 a positive in- 
teger. This relation introduces a nonlinearity in the de- 
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pendence of the mode frequencies n2, 713) = c\k\ 

on the parameters ni, 712 and 71,3, i.e. 

a;fc(ni, 712,713) _ ^ ^ nni 7rn2 27rrt3 ^ 



where 



fix) 



Li ' L2 ' cT 
-a;3 + ^/xlT^xfTxf) 



(3) 



As will be discussed in Sect. Ill M this is the main differ- 
ence between our approach and the free space calculation 
performed in 



II. THE COMMUNICATION RATE 

The communication rate R is the maximum number 
of bits per second that can be transmitted through the 
channel, and is given by the capacity (i.e. the maximum 
of the mutual information between the input and the out- 
put of the channel) divided by the transmission time T. 
The capacity can be estimated from a quantum mechan- 
ical analysis of the communication in which each symbol 
6, transmitted with probability density p(0), is associated 
with a quantum state p{6) of the Hilbert space H of the 
media used in the communication process (in our case the 
electromagnetic field). Without any constraints, the infi- 
nite dimensions of the input space Ti. of the system under 
consideration can accommodate an arbitrary amount of 
information and the capacity would diverge. Physically 
it is, thus, sensible to introduce an energy constraint on 
the accessible input states 0- In particular, we consider 
the following limit on the available average power (i.e., 
energy transmitted per unit time) 



P = 



Tr[gp] 

r 



(4) 



where p =^ J d9 p{9)p{9) is the average message sent 
through the channel (i.e. the electromagnetic field state 
at the input), and H is the Hamiltonian of the modes, 
i.e. 



k,e 



k 0,1, Mkf 



(5) 



fc,e 



with e =TE,TM and ak.c being the annihilation operator 
of the mode e with wave vector k and commutator 



Ofc, 



See' 5, 



kk' 



(6) 



In order to calculate the capacity of the channel under 
consideration, i.e. the maximization of the mutual infor- 
mation under the constraint (@J, we need the infinite di- 
mensional extension of the Holevo theorem @ which, 
in the noiseless case, gives an upper bound to the capac- 
ity in terms of the maximal input von Neumann entropy. 



As discussed in Ref. | 
so that the rate is 

R 



I, this upper bound is achievable. 



Sip) 
max — — 

p r 



(7) 



where S{p) = — Tr[plog2 p] is the Von Neumann entropy 
and the maximum is taken over all the possible density 
matrices p of the modes employed in the transmission, 
which satisfy Eq. Q). Notice that Eq. ((Tj) might be seen 
as an instance of the Holevo-Schumacher- Westmoreland 
theorem 7] , but for the case under consideration was first 
derived by Yuen-Ozawa ^2|. The maximization of Eq. Q 
under the constraint Q can be performed by means of a 
variational principle: the maximum is reached for the p 
that satisfies 



Sip) 
T 



where A and A' are the two Lagrange multipliers that 
derive from the power constraint ^ and from the nor- 
malization condition on p, respectively (the factor In 2 
has been inserted so that all calculations can be per- 
formed using natural logarithms). Using standard tech- 
niques (see for instance Q), it is possible to show that 
Eq. (|HJ) is satisfied by the density matrix 



-\oH 



ZiXo) ' 



(9) 



where Z{X) =Tr[e is the partition hmction and Aq 
is determined by the equation 



(10) 



Using this solution, the maximum rate in bits per unit 
time is finally given by 



R = 



1 



AnP 



lnZ(Ao 



(11) 



To obtain an explicit expression for R, we thus only need 
to evaluate the partition function Z(X) for the Hamilto- 
nian |SJ|. In the two following sections we will undertake 
such endeavor for two different communication scenarios. 



A. Multimode communication 

In this section we calculate the rate R when all the 
wave vectors that propagate in the positive longitudinal 
direction (from the sender to the receiver) are employed 
in the communication. 

Since modes with different k or different e are inde- 
pendent, the partition function factorizes in product of 
single mode partition functions Zk,ciX) so that 



InZ(A) =^lnZfc,,(A) , 



(12) 



fe,e 
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with 



with 



1 



-XhujL, 



(13) 



Substituting (|13|l into p2|) one can compute Z{X) by 
summing over the allowed values of ni, n2 and 713. Since 
we are interested in the stationary information rate, we 
should take the limit T — > 00 which allows the sum- 
mation over 713 to be replaced with an integral. Even 
with this simplification, the calculation is quite demand- 
ing and is postponed to the final paragraphs of this sub- 
section for the sake of readability. Here we consider the 
simpler high power/high cross-section regime, defined by 
the condition 
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AP 



> 1 



(14) 



where A = L1L2 is the cross sectional area of the waveg- 
uide. In this regime too the summations over ni and 712 
reduce to integrals and, apart from corrections of order 
1/7, Eq. ((T^ becomes 



InZ(A) 



27r3 



dx In 



1 



1 



-Xhcf{x) 



(15) 



where the volume integral must be performed on the sub- 
space V of positive Xj and f{x) is defined in Eq. (jSJ. The 
parameter g = 2 in Eq. H15|l counts the different species 
of modes, TE and TM in this case. It plays the same role 
as the polarization degeneracy in the free space propaga- 
tion of electromagnetic waves. By performing a change 
of integration variables and using the integral of Eq. (jA2|) 
of App. 13 Eq. reduces to 



InZ(A) 



gn 



AT 



240 h?\^c^ ■ 
Substituting this result in Eq. ifTIHl gives 

1 

A 



80 Pft3c2 



(16) 



(17) 



which, through Eq. (|11|) implies the following maximum 
rate 



R 



3 In 2 



gif^A 
80 c2 



(18) 



Numerical results. — When condition H14I) does not 
apply, the summations over rii and 112 in the definition 
of InZ(A) cannot be performed analytically. In this case 
we can resort to numerical evaluation of the rate. For 
the sake of simplicity we will consider a waveguide with 
square cross section, i.e. Li — L2. Remembering that 
not all the values of ni and 77.2 contribute both to the TE 
and to the TM modes, the summation of Eq. H12fl can be 
written as 



1 ry,^^ cT / -kXHc 



(19) 



where 

^ni ,712 (/^) 



dx In 



ni—l 



1 



1 — g-/3/(ni,7i2,x) 



(20) 



(21) 



Using Eq. I|ll|l , one can show that R vA/ c is a function 
only of the dimensionless parameter 7, defined in H14|l . 
In fact Eq. (|11|) becomes 



R 



ln2V^ 



W(/3o) 



(22) 



where (3q is the dimensionless quantity determined by the 
condition (|10|l . i.e., the solution of 



7 

TT 



(23) 



In Fig. 121 the numerical evaluation of the rate R is re- 
ported. Notice that in the limit 7 1, the solution H22|) 
approaches the asymptotic behavior 



Ry/A 



3 In 2 V 80 



9^^Y I 

' 'y 4 



(24) 



which corresponds to the high-power/high cross-sectional 
limit solution of Eq. (|18|l discussed previously. 

The dimensionless parameter 7 that identifies the on- 
set of the asymptotic regime for R^/A/ c has a relatively 
simple physical interpretation. From the density ma- 
trix, Eq. @, we see that the occupancy probabilities 
are highest for the lowest frequency (lowest photon en- 
ergy) modes. From the exact formulation for the parti- 
tion function, Eq. H12|l . we see that the triple-integral ap- 
proximation in Eq. H15|l will be valid when the occupancy 
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1000 
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FIG. 2: Numerical plots, a) Plot of the rate R given by 
Eq. 12211 as a function of the dimensionless parameter 7. b) 
Comparison between the same solution and the asymptotic 
behavior Rasym of Eq. 1241 : the ratio between these two quan- 
tities tends to one for 7 S> 1. 
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probabilities change very little between modes with adja- 
cent energy levels. The largest such photon-energy spac- 
ing occurs for the lowest-frequency modes, and, roughly 
speaking, is equal to hc/\/A, i.e., the photon-energy of 
the waveguide's cutoff frequency. If we concentrate all 
of the sender's average power into the lowest frequency 
mode, the resulting power spectral density will be ap- 
proximately P\/A/c, and hence 7 equals this spectrum 
measured in units of the photon energy hc/\/~A of the 
lowest-order mode. The condition 7^1 then guarantees 
the desired smooth behavior of the occupancy probabil- 
ities, because Eq. Q implies that high-photon-number 
occupancy of the lowest-order mode will, of necessity, be 
accompanied by similar occupancy of other low frequency 
modes. 

This analysis clarifies the regime of applicability of the 
approximation IjlHI) and underlines the importance of the 
quantity 7 in the definition of maximum rate. 



the wave vector k. We will specify this direction assign- 
ing the polar angle 9 — arccos(fc3/|fc|) and the azimuth 
angle — arctan(A:2/A:i). In terms of the discretization 
parameters rti, 712 and ns, these conditions become 



■mil 



n2 
ni 



L2 

— tan (p 



( lIHl 

V L2 



cos* 



(25) 



(26) 



B. Single directional mode communication 

In this section we calculate the rate R when the wide- 
band transmission is limited to using a single direction of 



where the nonlinear relation of Eq. ^ was used in de- 
riving Eq. (|26() . In this case, only those modes with k 
compatible with the chosen direction contribute to the 
partition function sum H12|l . i.e. 



00 00 



lnZ(A)= ^ J2 In 



713—0 ni ,712 — 1 



1 

I _ g-27rAn3/(T cos^ S) 



^n2,nitanipL2/Li ^ni ,2n3Li sin 0/[cT(l+tan2 cos^ 0] i 

(27) 



where, for the sake of simplicity only the TM mode has been considered and where the two Kronecker deltas take 
into account the conditions H25|) and H2t)|l . Again working in the high power/high cross-sectional regime (|14|l . the 
summations can be replaced with integrals and the Kronecker deltas become Dirac delta functions, so that we find 



InZ(A) ~ ^ / da; In 



1 



^ Xhcx^/ cos^ 



S {x2 — xi tan (f) S \ xi — X3 



smf 



y/l + tan^ if cos^ 9 



ttT cos^ 9 
12\h 



, (28) 



where V is again the subspace of positive components of 
Xj. Substituting this result into Eqs. (|10() and we 
now find that the maximum rate is 



R 




(29) 



If both the TE and TM modes of the chosen direction 
were used for the transmission, then a factor y/2 would 
appear in Eq. 

III. DISCUSSION 

In the previous section we calculated the maximum 
rate R for information transmission through an ideal 
metallic waveguide under an average input power con- 
straint. In particular we found that when the sender 



is using all the available modes, R scales as A^^'^P^^'^, 
as reported in Eq. H18|l . This scaling is reminiscent of 
the free space communication one 1]. The main differ- 
ence between the two cases is that, for waveguides all 
the positively propagating k vectors actually reach the 
receiver thanks to the reflecting walls of the waveguide. 
For frequency-w propagation over an L-m-long free space 
path between identical circular apertures of diameter D, 
there are approximately [D'^ ui / {8cL)]'^ low-loss propaga- 
tion modes, per polarization state In essence, the low- 
loss modes represent propagation angles that lie within 
the solid angle subtended by the receiver at the sender. 
These low-loss modes can be accounted for by introduc- 
ing a factor of V sin 9max {9max being the channel angular 
aperture as seen by the sender) in the free space rate — 
see Sect. VI A in 0. The lossy free space modes that 
do not satisfy the preceding angular subtense condition 
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can be used for communication, but their analysis re- 
quires inclusion of an accompanying noise source, which 
is mandated by the quantum theory of loss. As noted at 
the end of this section, finding the capacity of the lossy 
propagation channel is a considerably more difficult prob- 
lem. 

Apart from these physical considerations, a techni- 
cal difference between our calculation and the free space 
analysis is also evident. In deriving their result, the Au- 
thors of instead of using Eq. iQ, calculate the maxi- 
mum rate as 



R 



max 



E 

fc,e 



T 



(30) 



All the results discussed in this paper have been ob- 
tained in the lossless case, in which all the photons in- 
jected into the waveguide arrive to the receiver. In the 
presence of loss, the calculation procedure complicates 
noticeably: the capacity is no more simply given by the 
entropy of the initial state, but by the Holevo quantity, 
which is not known to be additive over successive uses of 
the channel ^ .IQj . Finally, we have not considered 
the presence of prior entanglement shared between the 
sender and the receiver. In this case the rate R can be 
doubled by using the super-dense coding protocol 
The calculation of the entanglement assisted capacity for 
the single directional mode case in the presence of loss 
was given in , following the procedure of . 



where 9 is the polar angle of the mode wave vector k 
and the maximum is performed over the mode states p^.e 
(e here counts the different polarizations of free space 
electromagnetic waves) . The presence of the term cos 9 is 
introduced in the sum to take into account the difference 
in the longitudinal speed of mode propagation. [In our 
calculation it is the nonlinear Eq. (0) that takes care of 
this]. Accordingly, the power constraint is calculated as 



P 



E 



■ COS( 



(31) 



IV. CONCLUSIONS 

In conclusion we have calculated the maximum com- 
munication rate for a perfect waveguide in the two 
regimes of multimode and single directional-mode com- 
munication. A comparison with the known results on free 
space communication schemes has been given. 

This work was funded by the ARDA, NRO, NSF, and 
by ARO under a MURI program. 



This procedure assumes implicitly that the maximum 
communication rate is achieved by a global state of the 
input modes which is unentangled over k. This assump- 
tion is correct as can be seen (at least for the waveguide 
communication protocol studied here) from the factor- 
ized form of the state in Eq. 1^. In order to compare the 
two approaches, we have calculated the maximum rate 
of the waveguide using Eqs. H3U|) and H31|) in place of Q 
and The results are, predictably, similar to the ones 
reported in Sec. Ill Al even though the numerical factor 
differs: in fact, for multimode communication and high 
power regime, we now find 



R = 



3 In 2 



120 c2 



(32) 



smaller than Eq. H18|l by a (3/2)^/'* factor, which derives 
from the particular choice of maximization of Eq. (|30() . 

If we consider the case of a single directional mode, on 
the other hand, R scales as cos 6ipi/2 (seeEq. this 
is the same result obtained in the free space propagation 
P, , apart from the cos 9 factor that takes into account 
the decrease in longitudinal propagation speed of the field 
due to the reflections at the waveguide walls. 



APPENDIX A 



In this appendix the integration needed for Eq. 1151) is 
given. 

After performing the change of variables yj = 2XhcXj 
{j — 1, 2, 3), the integral in Eq. p5(l becomes 



dyi / dy2 / dys In 

Jo JQ 



1 



— dx X 
2 70 Jo 



dy In 



1 - e-/(y) 
1 



(Al) 



where in the right hand term polar coordinates have been 
employed in the (2/1,1/2) plane. Changing to polar coor- 
dinates also in the plane spanned by (2a;,?/), the integral 
becomes 



(1 + sin ( 



dr r In 



1 



1 -e-^ 



120 



.(A2) 
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